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Abstract: We search for hidden mirror symmetries at large angular scales in the WMAP 7 year 
Internal Linear Combination map of CMB temperature anisotropics using global pixel based esti- 
mators introduced for this aim. Two different axes are found for which the CMB intensity pattern 
is anomalously symmetric (or anti-symmetric) under reflection with respect to orthogonal planes at 
the 99.84(99.96)% CL (confidence level), if compared to a result for an arbitrary axis in simulations 
without the symmetry. We have verified that our results are robust to the introduction of the galactic 
mask. The direction of such axes is close to the CMB kinematic dipole and nearly orthogonal to the 
ecliptic plane, respectively. If instead the real data are compared to those in simulations taken with 
respect to planes for which the maximal mirror symmetry is generated by chance, the confidence level 
decreases to 92.39(76.65)%. But when the effect in question translates into the anomalous alignment 
between normals to planes of maximal mirror (anti)-symmetry and these natural axes mentioned. We 
also introduce the representation of the above estimators in the harmonic domain, confirming the 
results obtained in the pixel one. The symmetry anomaly is shown to be almost entirely due to low 
multipoles, so it may have a cosmological and even primordial origin. Contrary, the anti-symmetry 
one is mainly due to intermediate multipoles that probably suggests its non-fundamental nature. We 
have demonstrated that these anomalies are not connected to the known issue of the low variance in 
WMAP observations and we have checked that axially symmetric parts of these anomalies are small, 
so that the axes are not the symmetry ones. 
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1 Introduction 

The simple and minimal ACDM standard cosmological model with the approximately flat (n s ~ 1) 
initial spectrum of adiabatic density perturbations predicted by the inflationary scenario of the early 
Universe is very successful in explaining the observed Cosmic Microwave Background (CMB) properties 
of temperature and polarization anisotropies pattern [1] , see also [2] and [3] for additional recent high-i? 
data from the Atacama Cosmology Telescope and the South Pole Telescope respectively, as mainly 
given by statistically isotropic and Gaussian fields. However, at the next level of accuracy several 
anomalies have been found using very different estimators. Though their contribution to the total rms 
values of CMB temperature anisotropy and polarization is small, they seem not to cope too easily 
with the statistical isotropy. At present it is not clear if these anomalies have a common origin. Also 
their statistical significance is under debate [4]. 

Namely, unlikely alignments of the quadrupole and the octupole were found in [5-9] for the WMAP 
first year release maps. Both quadrupole and octupole were shown to align with the CMB dipole [10, 11] 
(other unlikely alignments are described in [12-15] and a test for detecting foreground residuals in the 
alignment estimators is proposed in [16]). Also, it is found that the power in the CMB temperature 
anisotropy coming separately from the two hemispheres (defined by the ecliptic plane) is unlikely 
symmetric [17, 18]. It has been confirmed in the WMAP 3 year and 5 year release [19-21] and it is 
present in the COBE data as well, although with lower significance. The temperature power spectra of 
the opposing hemispheres are inconsistent at 3er to 4tr depending on the range of multipoles considered. 
The asymmetry has been detected in low resolution maps [17], both in angular and multipoles space, 
but it extends to much smaller angular scales in the multipole range S£ = [2,600] [20]. At large 
angular scales the hemispherical asymmetry has been tested for the first time in polarization maps in 
[22] making clear that this anomaly is mainly evident only in intensity at WMAP sensitivity. 

It has been suggested in [23] that an estimator built upon the point parity symmetry might be used 
as a practical tool for detecting foregrounds. In particular the authors consider whether the observed 
low CMB quadrupole in temperature could more generally signal odd point-parity, i.e. suppression of 
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even multipoles. However, they claim that WMAP 1st year data does not support parity preference 
beyond the meagre 95% confidence level. Later, it was found [24] that the parity symmetry in the 
temperature map of WMAP 3 and 5 year data is anomalous at the level of 4 out of 1000 in the range 
SI = [2, 18]. Other analysis in the WMAP 7 year data confirm the anomaly at same level for a slightly 
wider range Si = [2, 22] [25] and show that the anomaly has a characteristic scale around i ~ 20 [26] 
. See also [27, 28] for extensions of such analysis. 

Existence of these anomalies shows that CMB temperature fluctuations (as well as polarization 
ones) may contain some additional hidden information (pattern, "inscription") not evident from the 
standard analysis because it is masked by the dominant, statistically isotropic and Gaussian compo- 
nent. If related to low multipoles, this information can provide us with some completely unexpected 
knowledge about early stages of the evolution of our Universe including its very origin. However, it 
is not easy to extract it without having at least some crude idea of how it may look like. A hidden 
partial symmetry of CMB fluctuations, i.e. the existence of a subcomponent of CMB fluctuations for 
which this symmetry is an exact one, is one of the simplest variants of such hidden information. 

In particular, this symmetry can be the mirror one with respect to reflection in some plane. As 
was shown in Ref. [29], one physical mechanism leading to the appearance of a subcomponent of 
CMB angular temperature fluctuations possessing such symmetry is the existence of a non-trivial 
spatial topology: the non-compact T 1 topology, or the compact T 3 topology in which one of its 
topological scales is much less than two others (also called a slab space). For such spatial topology, 
the corresponding component of wave vectors of perturbations, say k z , is quantized, and then the 
subcomponent with the exact mirror symmetry with respect to reflection in the XY plane is produced 
by perturbations with k z — 0. However, it is evident that a similar effect can arise even in the absence 
of any non-trivial spatial topology but if, for some (may be pure statistical) reason, large-scale scalar 
(density) metric perturbations in the region of space including our last scattering surface have a 
significant component with a very weak dependence on one of the spatial coordinates. 

This type of partial symmetry was first searched for in Ref. [30] using the COBE data with the 
negative result placing a lower limit of the minimal present comoving size L of the torus of the order 
of 4 Gps. After that this limit has been raised much, and now is close to 2i?; s « 28 Gps for the cubic 
T 3 topology obtained using the most recent WMAP data and searching for identical circles in the sky 
[31] (Ris is the present comoving radius of our last scattering surface). Note, however, that a partial 
mirror symmetry can be seen in CMB fluctuations even for L > 2Ri s when there is no identical circles 
in CMB fluctuations. 

For this reason, in the present paper we have performed a new search for a hidden mirror reflection 
symmetry in the WMAP 7 year temperature maps developing estimators in the pixel domain and in 
the harmonic domain (the latter only for the full sky case). For generality, we look for a mirror 
reflection anti-symmetry, too, though it is less clear how such symmetry may arise. The paper is 
organized as follows. In Section 2 we define the used tools to perform the analysis reported in Section 
3. We discuss our results in connection with the issue of the low variance in WMAP observations in 
Section 4. Our conclusions are drawn in Section 5. Algebraic details for the equivalence between the 
estimators in harmonic and pixel domains are given in Appendix A. 
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2 Tools and estimators 



We consider the following pixel based estimators 

^5TV , ST 



Q(hi 
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where S^(hi) are defined in analogy to [30] as follows 
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where the sum is meant over the observed pixels, -/V p i x > $T/T(hj) are the temperature CMB anisotropics 
measured at the pixel pointed by the unit vector hj and fit is the opposite direction of hj with respect 
to a plane defined by hi, i.e. 



h k = hj - 2 (hi ■ hj)hi 



(2.3) 



We stress that only two out of the three global estimators given in Eqs. (2.1), (2. 2) are independent 
and that S ± is more sensitive to increased angular resolution than Q (see below). Eqs. (2.1) and (2.2) 
can be rewritten in the harmonic space as follows (see Appendix A for details): 
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where (£ + 1) C £ + = Yfm=-tPimW™\ 2 and £C e = X) m= _^ m K m | 2 with Ptm = 1 if + m ) is 
even/odd and pf m — if (£ + m) is odd/even. The estimators C e depend on the direction hi since the 
coefficients of spherical harmonics ai m are obtained in the frame in which the z— axis lies along the 
direction defined by hi. Notice that the quantity S + + S~ is always direction independent if computed 
over the whole sky, and it equals the total dispersion of the CMB fluctuations: 
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where Ct is the angular power spectrum of CMB anisotropies (see Appendix A). 1 

If different CMB anisotropy multipoles ai m are not correlated and statistically isotropic, then the 
average values of Q(hi), S + [hi) and S~(hi) are given by 



(Q(n t )) = -i-^Q 
i 

(S-(hi)) = 



(2.7) 
(2.8) 
(2.9) 



1 Note that our estimator Q differs significantly from the one adopted in the recent paper [32] which does not have a 
simple expression in the pixel space. 
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The r.h.s. of Eqs. (2.7), (2.8) and (2.9) are invariant under rotation, i.e. independent of the direction 
hi, as it should be. On the other hand, if CMB fluctuations contain a significant sub-component 
exactly symmetric (anti-symmetric) with respect to reflection in a plane with the normal vector rii , the 
estimator S~ (rij) (S + (ni)) will be anomalously small while the estimator S + (rii) (S~ (n*)) anomalously 
large. The maximum and minimum values for S~ are also maximum and minimum values for Q, 
respectively. 

3 The Analysis 

3.1 Data, Simulations and Results for the pixel based estimators 

We have considered the WMAP 7 yr Internal Linear Combination (ILC) map at Healpix 2 resolution 
-^Vgide — 16 (corresponding to maps composed by 3072 pixels, [33]), smoothed at angular resolution 
FWHM = 9.1285°. For each direction defined by the center of each pixel at such resolution, we have 
computed the estimators S + , S~ and Q defined by Eqs. (2.1) and (2.2). We have then repeated the 
same procedure for 10000 full sky CMB maps, extracted from the best fit of WMAP 7 yr ACDM 
model [1] at the same resolution and smoothed at the same angular scale as the WMAP map. Since 
random extractions are independent, CMB anisotropies are isotropically distributed on the sphere and 
the distributions of the estimators given in Eqs. (2.1) and (2.2) do not depend on the direction hi for 
the full sky case. 

Among all the directions connected to the HEALPIX pixelization at N^de = 16, for WMAP 7 
ILC map we find that the maximum value for S + corresponds to the minimum value for S~ , and vice 
versa. The minimum and maximum value of S + are plotted in Fig. 1 (dashed and solid vertical lines) 
on top of the distributions for the same quantities resulting from our 10000 Monte Carlo realizations 
(green histograms). The minimum values for S ± obtained by the WMAP 7 yr full sky ILC map are 
quite small compared to the probability distribution obtained by the Monte Carlo method, whereas 
their maximum values are not anomalous. The maximum and minimum values for Q do not exhibit 
anomalies with respect to the Monte Carlo distribution. We have explicitly checked that the probability 
distributions of S ± computed on WMAP 7 yr full sky ILC map over the 3072 directions associated to 
the HEALPIX N^c = 16 differ from the corresponding ones obtained from a random Gaussian CMB 
realization. The anomalous minima values for S are a consequence of the corresponding anomalous 
distributions obtained on WMAP 7 yr ILC map. 

It turns out that S + reaches its minimum value for the direction d\ defined by (0 — 107°, 4> — 264°), 
direction for which S~ is maximum. Instead, S~ is minimized along d 2 , defined by (9 = 42°, <j) = 260°), 
which also corresponds to the maximum of S + . Throughout the present paper, spherical coordinates 
are given in (9, 4>), being 9 and <f> the galactic co-latitude and longitude respectively 3 . For comparison, 
the direction corresponding to the hemispherical power asymmetry is (9 — 107°, <fi — 226°) [20], while 
that of the kinematic dipole measured by WMAP is (6 = 42° , (f> = 264° ) [34] that is very close to the 
direction d%- 

The probabilities related to the directions considered are summarized in Table 1. The most 
significant departure from symmetry is the S + , minimum along the direction d\ , followed by the S~ , 
minimum along the direction d%. The estimator Q shows no significant anomalies with respect to our 
Monte Carlo realizations. The secondary minimum in S + - S~ and Q - map is placed in the direction 

2 http:/ /healpix.jpl. nasa.gov 

^Spherical coordinates are sometimes provided in the literature in terms of galactic longitude and latitude, (i, b), 
related to the coordinates adopted in this paper by I = <f> and b = 90° — 9. The two directions d± and would therefore 
read (I = 264°, b = -17°) and (I = 260°, b = 48°) respectively. 
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Figure 1. Distribution of S + (left panel), S~ (middle panel) and Q (right panel) for the full sky analysis. 
For each estimator, vertical lines show the most anomalous values on the WMAP 7 yr ILC map: solid lines 
correspond to di, whereas dashed lines correspond to d\. The green histograms are obtained using 10000 
Monte Carlo realizations for S + (di) (left panel) and S~ (ds) (right panel). 




Figure 2. Distribution of S + (left panel) and S~ (right panel) for the masked sky analysis. S + is minimized 
along the same direction found over the full sky, di. The minimum of S~ points instead along direction d,3, 
with a small shift compared to the full sky result, d^. The blue histograms are obtained using 10000 Monte 
Carlo realizations for S + (di) (left panel) and S~(d,3) (right panel) with the galactic mask. 



characterized by (9 = 90° , <f) = 200°) - (9 = 171° ,(/> = 165°) -. The probabilities related to these 
second minima are 99.70% , 99.65% , 89.42% for S + , S~ , Q maps, respectively. 

Table 1. Probabilities (in percentage) to obtain a larger value for the considered directions and pixel based 
estimators. Numbers in brackets refer to directions obtained using the masked analysis. 



Estimator / Direction 


d x 


d 2 (d 3 ) 


S+ 


99.99 (>99.99) 


36.58 


s- 


18.78 


99.86 (99.99) 


Q 


2.35 


90.67 



We have also performed the same analysis by masking the galactic plane and considering in the 
sums only the contribution from those pixels which are not masked together with its opposite. The 
factor l/A^pix in Eqs. (2.1,2.2) takes into account only these pixels. The comparison of the anomalous 
directions between full sky and masked sky shows that they are stable at the chosen resolution (see 
Fig. 3), as only a minor deviation appears for the minimum of S~ along the direction d 3 of coordinates 
(6 = 39°, 4> = 260°). However, the level of anomaly increases for the masked sky for both S + and S~ 
(see Table 1 and figure 2). 
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Figure 3. Sky projection of the relevant directions: light blue for di; red for d,2\ blue for ^3; yellow for the 
N/S asymmetry; green for the kinematic dipole 

So far we have reported probabilities based on a comparison between real data and simulations 
without symmetries with respect to a fixed axis. If instead the real data are compared to the mini- 
mum of S ± across the sky in each of 10000 random simulations, the statistical significance drops to 
92.39(76.65)% for S+ (S~) on the WMAP 7 yr ILC map, as can be seen in Fig. 4; 14.77% of the 
simulations lie to the left of the the minimum value for Q on the WMAP 7 yr ILC. On the basis of this 
second comparison of absolute minimum values of estimators on real data vs on simulations across the 
sky for each random realization, the anomaly translates into the question of an anomalous alignment 
between the normals to the planes of maximal mirror symmetry of real data with known directions 
such as the normal to the ecliptic plane and the cosmological dipole. 

In Fig. 5 we provide now the full sky score map for Q, i.e. Q{hi) in Eq. (2.1). Very similar maps 
are the full sky score maps for S ± , i.e. S (hi) in Eq. (2.2) (where the two anomalous directions are 
shown in grey), which are given by ±Q shifted by the value of the variance, which is represented by 
a monopole in these maps since it is direction independent. The full sky maps are even under (point) 
parity symmetry. The score maps are shown in Fig. 5 also for the masked case, although the masked 
score maps for 5' ± and Q are all different. 
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Figure 4. Distribution of S + (left panel), S~ (middle panel) and Q (right panel) for the full sky analysis. 
The most anomalous values on the WMAP ILC map (vertical line) is compared to two sets of MC simulations: 
in violet the estimators for a fixed axis (as in Fig. 1) and in pink the minimum value across the sky computed 
independently for each sky. 



S+ full sky 3+ masked sky 




Figure 5. Map of S + (top), S (middle) and Q (bottom) for the full (left) and masked (right) case. 
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3.2 Data, Simulations and Results for the harmonic based estimators 

In order to study the angular scale dependence of the two anomalous directions found through the 
application of the pixel based estimators, we have considered the harmonic version of the same esti- 
mators, given in Eqs. (2.4) and (2.5). 

We have smoothed the WMAP 7 yr full resolution ILC map at angular resolution FWHM = 
2.2903° and reconstructed at lower resolution N si d c = 64. We have then extracted 10000 CMB 
realizations at the same resolution iV S idc = 64 from the best fit of the WMAP 7 year ACDM model 
and smoothed at the same angular scale as the WMAP map. For each of the two directions found 
in the pixel-based analysis, we have computed S + , S~ and Q as in Eqs. (2.4) and (2.5), for both 
WMAP and simulations, substituting the indefinite sum over a multipoles with a definite one between 
£ min = 2 to £ max = 80. By comparing those, we have computed the probability of WMAP with respect 
to the random and independent CMB Gaussian extractions. Results are plotted in Fig. 6. The most 
anomalous case turns out to be the S + estimator for the direction d± and for € max > 15 where the 
probability of getting such a low value is 0.01%, see upper right panel of Fig. 6. The estimator S~ 
for the direction di is, on average over I, anomalous at the level of ~ 0.3% for i max > 4 with a peak 
reaching 0.13% for £ mgcK = 23 (see middle left panel of Fig. 6). 

It is important to assess the weight of the low multipoles in quantifying the anomaly of the two 
directions. To do so, we have analyzed various cases not including the first multipoles from the analysis. 
More precisely, we consider three different cases, by excluding from the analysis the quadrupole, both 
quadrupole and octupole and the first 9 multipoles (the sum starts from £ = 11 in the latter case). 
The corresponding results are shown in Fig. 7. For the d^ (symmetry) direction, the quadrupole 
has an impact in reducing the anomaly and by cutting the first 9 multipoles the anomaly almost 
disappears. Thus, it seems that the quadrupole makes some contribution to the effect of anomalous 
mirror symmetry, but it is not the dominant one. The reported anomaly here is therefore not just 
given by the anomalous alignment of the quadrupole and octupole found in [10]. On the other hand, 
for the g?i (anti-symmetry) direction, there is a noticeable change when cutting the low multipoles, 
but an anomalous behaviour still persists. 

On summarizing the first three panels of Fig. 7, it is clear that the anomaly for S + has a 
contribution from intermediate multipoles, 10 < £ < 35, whereas S~ anomaly is mainly due to low 
multipoles and almost disappears for I > 10. 

Since the direction of our d 2 axis is close to that of the CMB dipole which is mainly kinematic 
and local (compared to cosmological scales), the question arises naturally if it may be caused by 
an imperfect dipole subtraction or, more generally, if it is somehow related to the dipole which is, 
in particular, axially symmetric. Also it is interesting to check if there exists some traces of axial 
symmetry with respect to the axes we found (and which define planes of partial mirror symmetry 
orthogonal to them). For these purposes, we repeat our investigations excluding now the axially 
symmetric component with m = from the estimators S + and S~ using their multipole representation 
(2.5) with respect to the di axis, as well as to the d\ axis, too (though in the latter case there is no 
reason to expect a significant amount of axial symmetry). This method is similar to that used in 
the recent paper [35] in the case of the partial parity symmetry. Results are presented on Fig. 5, 
too. It is clearly seen that previous curves for these estimators are practically unchanged (in fact, the 
anomaly even becomes slightly more for the di direction). This clearly shows that axially symmetric 
contributions to the anomalies found are small and that the directions d\ and di may not serve as 
symmetry axes. 
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Figure 6. Percentage versus the multipole number ^ ma x- For each ^ max and for each estimator, we plot the 
percentage to obtain a larger value w.r.t. the WMAP ILC one. Left column is for the d,2 (42°, 260°) direction 
and right column is for the di (107°, 264°) direction. First row for S + , second row for S~ and third row for 
Q. 

4 Discussion 

One might argue that the weight of the variance term (direction independent) in our estimators S ± , Q 
in Eqs. (2.1,2.2) connects our reported anomalies on the mirror-parity directions to the known problem 
of the low variance in WMAP maps. We have computed the low variance of WMAP 7 yr ILC at (with 
a Gaussian smoothing of 9.1285 degree) against MC simulations based on the WMAP 7 yr ACDM 
best-fit in Fig. 8: although we cannot state that the WMAP variance is anomalous at a significant 
statistical confidence level, it is indeed smaller that the average value of the MC distribution, as shown 
in Fig. 8. 
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In order to try to disantangle the issue of the low WMAP variance from our findings, we also 
compare our anomalous values with 10000 MC simulations based on the WMAP 7 yr temperature 
power spectrum (publicly available at the NASA Lambda web site) which contains the small values for 
the quadrupole and octupole, rather than the theoretical ACDM theoretical best-fit. The comparison 
of our results against this second MC shows that the reported directions for are again anomalous, 
respectively at 99.50% and 99.17% confidence level, but a little less than with respect to the MC 
based on the WMAP 7 yr theoretical best-fit, in agreement with the WMAP low variance shown in 
Fig. 8. The most intriguing result is that the minimum and maximum values of Q become much 
more anomalous with respect to this second MC, respectively at 99.83% and 99.24% confidence level: 
this fact gives an alternative perspective to our results, now connecting our reported anomalies to the 
correlation across the plane, as shown in figure 9. Table II displays the probabilities with respect to 
the second MC based on the WMAP 7 yr temperature power spectrum. 

Table 2. Probabilities (in percentage) to obtain a larger value for the considered directions and pixel based 
estimators with respect to the MC with 10000 realizations of the WMAP 7 yr temperature power spectrum 



Estimator / Direction 


di 


d 2 


S+ 


99.50 


6.15 


s- 


2.47 


99.17 


Q 


0.17 


99.24 



We have corroborated our analysis in pixel space by the analysis in the harmonic domain: we 
therefore present plots analogous to Fig. 4, now with respect to the temperature WMAP 7 yr power 
spectrum, in Fig. 10. The plots confirm the analysis in the pixel domain, i.e. that our reported 
anomalies presented in section 2 are not related to the low variance problem of WMAP 7 yr map and 
that the correlation across d\ and d 2 are anomalous with respect to this second MC. This analysis 
clarifies the anomaly reported with respect to MC simulations which keeps fixed the axis of (anti-) 
symmetry to the anomalous ones. We do not repeat this analysis with respect to simulations which 
vary the direction in which the maximal (anti-) symmetry is generated. 

5 Conclusions 

We have searched for hidden mirror symmetries in the WMAP 7 year ILC temperature map. For this 
aim, we have developed new global pixel estimators and their relative expressions in the harmonic 
domain. We have found that two different axes exist for which the CMB temperature anisotropy 
pattern is anomalously symmetric (anti-symmetric) at the 99.84 (99.96)% confidence level, if compared 
to a result for an arbitrary axis in simulations without the symmetry. The angle between them is about 
65°. These axes correspond to the minimum and maximum of CMB temperature correlations with 
respect to reflection in a plane. We have verified that the above mentioned probability of the anomalies 
and the position of the axes are robust to the introduction of the galactic mask. The axis d\ for which 
the amount of mirror anti-symmetry is maximal is nearly orthogonal to the ecliptic plane, while the 
axis d 2 with the maximal amount of mirror symmetry is close to the kinematic dipole direction. It is 
interesting to note that the d\ axis is also close to the one which identifies the N/S asymmetry [21]. 
Removal of CMB temperature multipoles with m = with respect to these axes does not change the 
anomalous behaviour of our mirror symmetry (anti-symmetry) estimators showing that the symmetric 
(anti-symmetric) patterns found are not axially symmetric with respect to the axes d± and d 2 . Our 
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findings differ from those previously published on mirror symmetries [36, 37] using different estimators. 
Neither of those axes is close to the preferred direction found in the recent paper [38] for a model of 
anisotropic dark energy. 

If instead the real data are compared to the minimum of S ± , Q across the sky in each of 10000 
random simulations, the statistical significance drops to 92.39(76.65)% for S + (S~) on the WMAP 7 
yr ILC map, as can be seen in Fig. 4; 14.77% of the simulations lie to the left of the the minimum 
value for Q on the WMAP 7 yr ILC. On the basis of this second comparison of absolute minimum 
values of estimators on real data and simulations, the anomaly we find translates into the question of 
the anomalous alignment between the normals to the planes of maximal mirror symmetry of real data 
with known directions such as the normal to the ecliptic plane and the cosmological dipole. 

To clarify the nature of these anomalies, we have investigated the weight of the low multipoles 
in these new anomalies by exploiting the estimators in multipole space: by removing the quadrupole 
and octupole from the analysis, we have shown that the (anti-) mirror symmetries reported here are 
novel results and are not simply a different perspective of the previously discovered anomalies in the 
quadrupole and octupole in intensity power spectrum and multipole vectors [10]. By cutting I < 10, 
the statistical significance of the anomaly in both directions drastically drops. However, whereas 
the anomaly in S + (its anomalously low value of for the d\ direction) remains though becomes less 
pronounced, the anomaly in S~ for the di direction almost disappears. In our opinion, this different 
behaviour might be not completely surprising: the direction d\ is close to that of the N/S asymmetry 
which has been shown to be present also at intermediate multipoles of the order of I ~ 100 [21]. This 
suggests that the (anti-symmetric with respect the mirror reflection) anomaly in S + and the N/S 
asymmetry may have the same and probably secondary (even local) origin, that does not prevent it to 
be a real physical, non-systematic effects. On the opposite, the symmetry anomaly which we find in 
S~ is due only to low multipoles and therefore it may have a cosmological or even primordial origin. 
Its anomalous alignment with the cosmic dipole may be a sign of their common origin, too. Whatever 
the origin of these anomalies is, these two new anomalies differ in this way. 

We have shown that these hidden mirror symmetry properties have little to do with the low 
variance in WMAP maps. To gauge away the issue of the low WMAP variance from our results, we 
have also compared our findings with 10000 MC simulations based on the WMAP 7 yr temperature 
power spectrum, which contain the low quadrupole value. With respect to this second MC, the 
anomalies for do not disappear, whereas the anomaly for Q becomes more evident and connects 
our reported anomalies to the correlation across the plane. 

While the present paper was being prepared for submission, the article [32] appeared where a 
search for the mirror (anti-)symmetry in CMB fluctuations was also performed. In contrast to our 
paper, the analysis in [32] is restricted only to low multipoles up to i max = 7 for an estimator defined in 
the harmonic domain which is similar to Q, but with a I dependent weight. The statistically significant 
anomaly reported in [32] is an anti-symmetry anomaly in the direction (0 = 109°, (f> = 266°) which is 
close to our direction d\. We stress that we do not report any anomaly for Q in the WMAP 7 yr ILC 
map, but the different estimator and different data analysis adopted here and in [32] do not allow a 
direct comparison of the results. 

So, summarizing, we have found two planes, with respect to reflection in which CMB angular 
temperature fluctuations possess the most pronounced and anomalous amount of partial mirror sym- 
metry or anti-symmetry. From these two cases, the mirror symmetry one for which the corresponding 
plane has the normal direction di is more promising for the purpose of extracting new cosmological 
information about the present and early Universe since it is almost entirely due to low multipoles I. 
It is clear that the amount of existing data is not sufficient to make more definite conclusions about 
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the origin of this anomaly. The forthcoming release of WMAP 9 year data and then the independent 
observations from Planck will shed more light on these new mirror anomalies. Additional independent 
observations or different data analysis will be cross-checks of the importance of possible systematics, 
as beam asymmetries, foreground residuals and subtleties related to the kinematic dipole in the mirror 
symmetries found in this paper. 
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A Multipole expansion of the pixel based estimators 

In this Appendix we give the relevant algebra to Eqs. (2.4-8). We consider the expansion of temper- 
ature anisotropics in spherical harmonics 




(A.l) 



with the following convention 




(A.2) 



The convention on the averages is N pi ^ Ylj=T = ^ I By using the multipole expansion, Eq. (2.2) 
and the completion property of the spherical harmonics in Eq. (2.2), we obtain: 




(A.3) 



(A.4) 



where we have used: 




(A.5) 
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The last term in Eq. (A. 4) is simply — Q and can be rewritten as: 

J dSl^(0,<p)^(«-0,<p) = ^^a tm ai, m ,(-lf^\ m (0,<p)Y; lmt (e,<p) (A.6) 

Ira I'm' 

= £(-i)' +m kJ 2 (A.7) 

im 

= ^2(t + l)Ct-iCT, (A.8) 

t 

where Cf , CJ are defined below Eq. (2.5). 
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Figure 7. Percentage versus the multipole number £ max excluding 1 = 1 (first row), £ = 2,3 (second row), 
I < 10 (third row) and excluding the m = modes (fourth row). For each £ max and for each estimator, we 
plot the percentage to obtain a larger value w.r.t. Left column is for S + in its anomalous direction di, right 
column is for S~ in its anomalous direction of d 2 . 
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Figure 8. Value of the variance of the WMAP 7yr ILC map (thick error bar) compared with its MC 
distribution based on 10000 realizations of the WMAP 7 yr ACDM best-fit. 
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Figure 9. Value of the estimators S , Q in the WMAP 7yr ILC map (thick error bars) compared with its 
MC distribution based on 10000 realizations of the WMAP 7 yr ACDM best-fit (green histograms) and of the 
WMAP 7 yr temperature power spectrum (purple histograms). 
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Figure 10. Percentage versus the multipole number £ max with respect to the WMAP 7 ys observed power 
spectrum. For each ^ ma x and for each estimator, we plot the percentage to obtain a larger value w.r.t. the 
WMAP ILC one. Left column is for the d 2 (42°, 260°) direction and right column is for the di (107°, 264°) 
direction. First row for S + , second row for S~ and third row for Q. 
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